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0.1 Introduction 


This thesis is devoted to several problems posed for special classes of matrices (such as GKK) 
and solved using structured matrices (such as Toeplitz) belonging to that class. 

The topic of Chapter 1 is GKK r-matrices. This notion was introduced in the 1970’s 
as a response to the Taussky unification problem posed in the late 50’s, which is discussed 
in Section 1.1. Section 1.2 is devoted to four conjectures proposed in the 1970’s-1990’s 
on the stability of GKK r-matrices. They are all disproved in Section 1.3 using GKK r- 
matrices with additional structure (Toeplitz and Hessenberg). Further properties of the 
counterexample matrices, which are themselves not important for disproving the stability 
conjectures but seem to be worth analyzing, are taken up in Section 1.4. Section 1.5 contains 
a brief discussion of open problems related to the topics of earlier sections. 

Chapter 2 is centered around the following problem: given a collection of matrices (A a ) 
in a special class (such as totally nonnegative) bounded in some matrix norm and such that 
the spectrum of A a lies outside a disk of fixed radius with center at zero, determine whether 
the collection (A” 1 ) is bounded in the same matrix norm. For any matrix norm, the answer 
is yes for matrices of bounded order, as is shown in Section 2.1. The next sections all deal 
with collections of matrices of unbounded order and the ‘simplest’ oo-norm, the choice also 
motivated by applications. It is shown in Section 2.2 that the answer is still yes for totally 
nonnegative Hermitian matrices. However, the answer is no for positive definite Hermitian 
matrices. Section 2.3 contains a pertinent counterexample and a variation of it both based 
on the Hilbert matrix. A counterexample for the class of Hermitian Toeplitz matrices is 
obtained in Section 2.4. Finally, an interesting question of the same type arising in spline 
theory is discussed in Section 2.5. 

0.2 Notation 

The following conventions are used throughout the thesis. 

To make a clear distinction between equality and equality by definition, the latter is 
denoted by : = . The symbol # denotes the cardinality of a set. The set {1,..., n} for n G IN 
is denoted by (n). For p, q e 7Z, let 

if p < q 
otherwise 


p:q := l {p,p+l,...,q- l,q} 


For any x G IR, set 


rr + : = 


x if x > 0 

0 otherwise 
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The linear space of (column-)vectors with n entries in (D is denoted by <D n . The jth vector 
of the standard basis of <D n , i.e., the vector with 1 at the jth position and zeros elsewhere, 
will be denoted by e r The linear space of all n x n-rnatrices with entries in C is denoted 
by C nxn . If the order of a matrix is not clear from the context, it will be indicated by the 
subscript n x n or simply n. The symbol / stands for the identity matrix of appropriate 
order. Both the zero matrix and the number zero are denoted by 0. 


0.3 Basic matrix notions 


Given a matrix A e C nxn , i e t A(a, /3) denote the submatrix of A whose rows are indexed by a 
and columns by 0 (a, 0 G (n)) and let A[a, 0] denote det A(a, 0) if = #/3. For simplicity, 
A(a) will stand for A(a, a) and A[a] for A[a, a]. By definition, A[0] := 1. Elements of A are 
denoted by a(i,j). A block diagonal matrix with (square) blocks A and B will be denoted 
by 

diag(A,5)(:= ^ q ^). 

The spectrum of A, i.e., the multiset of its eigenvalues (with each eigenvalue repeated 
according to its multiplicity), is denoted by cr(A). The spectral radius of A is denoted by 
p(A)(: = max |er(A)|). 

The inequality A > 0 (> 0) means that A is entrywise nonnegative (positive). A > B 
(. A > B ) means, by definition, that A — B > 0 (A — B > 0). 

A norm [| • || on the space (C nxn is a matrix norm if it satisfies the inequality 

||AB|| < ||A|| • ||5|| \/A,Be€ nxn . 


A matrix norm j| • || 0 is the operator norm subordinate to the norm || ■ || on C n if 


|| A|| = sup V A G C nxn 

^ec n \{o} |p ; || 


In particular, the p-norm 
the p-norm 


v 



(1 < p < oo) on C nxn is the operator norm subordinate to 


(E”=i Hm 1 '" 

maxj=i ) ... )Tl |u(i)| 


if 1 < p < oo 
if p — oo 


on the space C n . The condition number of an invertible matrix A (for the norm || • ||) is the 
product ||A|| ■ HA^ 1 !!. 

A matrix T e C nxn is Toeplitz if it has the form T —:(r(i — j))fj =i for some (r(i))”T 1 1 _ n . 
A matrix is Hessenberg if the entries on its first subdiagonal are all equal to 1 and the entries 
below that subdiagonal are zero. 
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Chapter 1 

Eigenvalues of GKK matrices 


1.1 Taussky unification problem 


A matrix A is called totally nonnegative if A [a,/?] > 0 for all a, (3 G (n) with j^a = #/?. A 
is called an M -matrix if A = rl — P where P > 0 and r > g(P). For more than a dozen 
other ways to define M- matrices, see |2J. 

A matrix A is called a P-matrix if A [a] >0 for all a C (n). A is said to be sign-symmetric 

if 


A [a, /3]A[/3, a] > 0 \/a,/3 e (n), = #f3. 


A is called weakly sign-symmetric if 


A [a, 0\A\p, a] > 0 

for all Wa,/3 G (n) , with 

= # = (1A) 

The minors A [a,/?] with the property (II .1 1) are sometimes called almost principal. 

Weakly sign-symmetric P-matrices are also called GKK after Gantmacher, Krein, and 
Kotelyansky. 

Let 

1(A )miner (A) n IR, 

with the understanding that, in this setting, min0 = oo. A matrix A is called an u-matrix 
if it has eigenvalue monotonicity in the sense that 

l(A(a,a)) < l(A(/3,/3 )) < oo whenever (n). 

A is a r-matrix if, in addition, 1(A) > 0. 

Hermitian positive definite, nonsingular totally nonnegative, and M-matrices all en¬ 
joy positivity of principal minors, weak sign symmetry, and eigenvalue monotonicity. In 
fact, these properties were singled out as a response to the ‘unification problem’ for the 
above-mentioned three classes of matrices that stems from a research problem posed by 
O. Taussky [223 • 1 

1 0. Taussky pointed out in mi that similar theorems were known for some positive matrices and for positive definite 
Hermitian matrices, for which the then available proofs were different, and asked for a unified treatment of both cases. She gave 
four examples of such similar theorems, two of which illustrate common properties of totally nonnegative and positive definite 
Hermitian matrices. The term ‘Taussky unification problem’ was later taken to mean a much wider class of problems. 
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The fact that Hermitian positive definite matrices are sign-symmetric P-matrices (a prop¬ 
erty stronger than being GKK matrices) is standard. The eigenvalue interlacing property of 
Hermitian matrices (sa or, e.g., 0 p.59]) implies their eigenvalue monotonicity. Directly 
from the definition, nonsingular totally nonnegative matrices are sign-symmetric with non¬ 
negative principal minors. Their eigenvalue monotonicity was proved by Friedland |TI|j. This 
fact and the spectral theory of totally nonnegative matrices show that all principal minors 
of a nonsingular totally nonnegative matrix are in fact positive. 

The Perron-Frobenius spectral theory of nonnegative matrices shows that M-matrices are 
P- and uj- matrices. The weak sign symmetry of M -matrices was proved by Carlson [8]. 

1.2 Stability conjectures on GKK r- matrices 

Yet another property shared by Hermitian positive definite, totally nonnegative, and ill- 
matrices is their positive stability. To recall, a matrix is called positive (negative) stable if 
its spectrum lies entirely in the open right (left) half plane. In the sequel, the term ‘positive 
stable’ will be usually shortened to simply ‘ stable’. 

Hermitian positive definite and totally nonnegative matrices are obviously stable (hav¬ 
ing only positive eigenvalues), while the stability of M-matrices follows from the Perron- 
Frobenius theory. 

The natural question arising from this observation was whether some combination of the 
properties from Section 1.1, viz., positivity of principal minors, weak sign-symmetry, and 
eigenvalue monotonicity, is sufficient to guarantee stability. (None of those properties alone 
is sufficient, which can be checked by simple examples.) 

Carlson |HJ conjectured that the GKK matrices are stable and showed his conjecture to 
be true for n < 4. 

Engel and Schneider PS] asked if nonsingular r-matrices or, equivalently, cu-matrices all 
whose principal minors are positive (see Remark 3.7 in im are positive stable. Varga [221 
conjectured even more than stability, viz. 

I arg(A - l(A))\ < - - - VA G a{A). 

2 n 

This inequality was proven for n < 3 by Varga (unpublished) as well as by Hershkowitz and 
Berman 12 a and for n = 4 by Mehrmann eg. In his survey paper ea, Hershkowitz posed 
the weaker conjecture that r-matrices that are also GKK are stable. 

The above conjectures were plausible not only because they were verified for matrices of 
small order, but also due to the following two theorems. 

The first indicates that there is a certain ‘forbidden wedge’ around the negative real axis 
where eigenvalues of a P-matrix cannot lie. (The angle of the wedge depends on the order 
of the matrix.) 

Theorem (Kellogg [I2B1J). Let A G (D nxn be a P-matrix. Then 

I arg(A)| <ti -VA G <j(A). 

n 
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The second shows that sign symmetry together with positivity of principal minors is 
sufficient for stability. 

Theorem (Carlson ^|). Sign-symmetric P-matrices are stable. 

Carlson’s elegant proof employs the Cauchy-Binet formula to show that A 2 is a P-matrix 
whenever A is sign-symmetric, diagonal stable scaling property of P-matrices to conclude 
that DA is stable for some diagonal matrix D with positive diagonal, so that the homotopy 
S(t) :=((1 — t)D + tI)A preserves sign symmetry as well as positivity of principal minors, 
hence the eigenvalues of S(t) cannot cross the imaginary axis as t runs from 0 to 1. 

1.3 Counterexample 

As is shown in |2H] . none of the proposed conjectures is true. Described below is a class of 
GKK T-matrices which are not even nonnegative stable, i.e., they do have eigenvalues with 
negative real part. This class consists of Toeplitz Hessenberg matrices A n ^ yt of order n that 
depend on two more parameters k G IN and t G 1R. In what follows, it will be shown that 
A n j %t is a GKK r-matrix for any t G (0,1) and that A 2 k+ 2 ,k,t is unstable for sufficiently large 
k and sufficiently small positive t. 

Let A be an infinite Toeplitz Hessenberg matrix with first row (a 0 ,ai,...) and let d n 
denote its leading principal minor of order n. By the Laplace expansion by minors of the 
first row, 

n 

d n = ^(-l) n_1 a i _id n _ i , n G IN. (1.2) 

3 = 1 

This is, in effect, an invertible lower triangular system for the a 3 . So, for an arbitrary 
sequence (d\, d 2 , ...), there exists exactly one Toeplitz Hessenberg matrix having these as its 
leading principal minors. 

With this, let A^^t be the Toeplitz Hessenberg matrix with leading principal minors 

d n = t (.n-k-i)+, nG ] n, 

for some t G (0,1) and k G IN. Then equation (11,21) becomes 

n —1 

J^(-l ) j a/ n - j - k ~ 2) + = t {n - k - 1)+ , n G IN. (1.3) 

3=0 

Let A n> k,t be the leading principal submatrix of A 0 of order n. 

Two observations are immediate: The first k + 1 entries in the sequence (di, d 2 ,...) equal 
1 and dk +2 = t, hence a 0 = 1, a 3 = 0 for 1 < j < k, and (ik+\ = (—l) fc (l — t). Secondly, the 

matrix A^^t is Toeplitz, so all its principal minors indexed by j consecutive integers equal 
fU-k- 1)+_ 
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1.3.1 On two characteristics of A^^t 

As is well known, associated to any infinite Toeplitz matrix T — j))fj =0 is its symbol 

S{T), i.e., the Laurent series 

OO 

S{T)(s):= 

j = - oo 

It is also convenient to introduce the Taylor series T>{T) 


V(T)(s ,\) :=Y,Dj( A)** 

3=0 


where 

n ( X ) — / 1 if 3 = 0 

— \ det(T(l:j) - XIj) if j G IN ' 

To avoid cumbersome notation, Dj( A) will be used in the sequel, but one should keep in 
mind that Dj( A) also depends T, i.e., for T = on the parameters k and t, so should 

have been denoted by Dj( A, k,t). 

Proposition 1. 


5(Ao 0 ,k,t)(s) 

X>(Aoo,fc,t)(s, A) 


1 + ts 

1 + (1 - t) Efii si 

l+s(A-i) + A(l 


(1.4) 

(1.5) 


Proof. From & 

i 

ai = (-i ytv-v++ (1.6) 

3 =1 

for all l G IN. Let <E>(s) := <S(Aoo,fc,t)(s) — 1/s. Replacing each a v j G IN, by its expansion (II.hi) 
and collecting terms with the factor tb- fc )+ f or each j, one obtains 

OO OO 

$(s) =ao- ®(s)Y t{j ~ k ~ 1)+ h s ) j + Y tU ~ k)+ (~ s y- 

3 = 1 i =1 


Recall that ao = 1 and compute the sums in the right-hand side. This yields 

1 + ts — (1 — £)(— s) k+2 1 + ts — (1 — t)(s) k+1 

(1 + ,)(1 + U ) = (! + ,)(! +to) ’ 


hence 

S(A n kt)(s) = $(s) H— 

s 

which, after canceling (1 + s), gives (11.41) . 


(1 + s) (1 + ts) 
s(l + ts — (1 — t)(—s) k+2 ) ’ 


(1.7) 
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Now observe that, similarly to dm expansion of the determinant D n ( A) by its first row 
gives 

TO—1 

(do — \)D n _i+ '^^(—iyajDj(\) = D n (\), n 6 IN. (1-8) 

j =i 

That is, the coefficient of s J in X>(A»,fc,t)( s > A) equals the coefficient of s J ' -1 in the product 
of I?(^4oo ,k,t)(s,X) and $(— s) — A whenever j > 1. With the zeroth coefficient taken into 
account, this means 

t k,t){s, A) = s(<f ) (—s) — A) 'D{A 00 j t f)(s, x ) + 1) 

which, together with (11.71b implies (11.511 . □ 

Corollary. The polynomials Dj satisfy the recurrence relation 

k +2 

Dj (A) + (A - t)Dj_ i(A) + A(1 ~t)J2 D J-i ( X ) = 0 Vj > k + 2. (1.9) 

1=2 

Proof. Compare the coefficient of s- 7 for j > k + 2 in the numerator (zero) with that in 
the product of r D{A 00 ^,t) and the denominator in (11.51) . □ 

These results will be revisited in subsection 1.3.3. Let us now show that A^^y is GKK. 

1.3.2 A n k t are GKK 

Since A n ^ it is Hessenberg, the submatrix A n> kj((n)\i:i+j—l) is block upper triangular if 
l<i<i+j — 1 < n, so 

A U) k,t[ a u P\ = A n:k j[a}A n:ktt [/3] whenever i < j - 1 for all i e a, j G ( 3 . (1.10) 

This shows that A n ^ it is a P-matrix. 

Fortunately, there is no need to verify the weak sign symmetry of a P-matrix by computing 
its almost principal minors. Instead, one makes use of the following remarkable fact. 
Theorem (Gantmacher, Krein jlYt p.55], and Carlson |S|) . A P-matrix A is GKK 
iff it satisfies the generalized Hadamard-Fisher inequality 

A[a]A\f3\ > A[a G f3\A[a V\ 0\ Va,/3C(n). ( 1 . 11 ) 


Since 0 < t < 1 and 

(x + y) + + (x + z) + < x + + (x + y + z) + Vx Vy,z> 0, 

one obtains 

A n ^k,t\i-i~\-j 1 ] • A nk f\l.l-\-m 1 ] 

— f;U-k-l) + +(m-k-l)+ > ^{l+m-i-k-l) + +(i+j-l-k-l) + •£ ^ < i+j- 

A n) k,t\i-Ipim 1 ] ’ A n k ff ,i~\-j 1 


( 1 . 12 ) 


Together with cm cm shows that A n ^ satishes cm if a, P are sets of consecutive 
integers. 

To prove cm in general, first make a definition. Call the subsets a, P C (n) separated 
if \p — q\ > 1 Wp G a, q G [3. Suppose ct, Pi, ..., f3j C (n) are sets of consecutive integers, Pi 
(i — 1 ,..., j) are pairwise separated, and 

for any i — 1,..., j, there exist p G Pi and q G a such that \p — q\ < 1. (1.13) 

Then A nik>t , a, and P : = U { =1 Pi satisfy (11.111) . Indeed. (11.111) holds for a and Pi. If 1 < l < j, 
then, assuming cm for a and 7/ := U* =1 /?i, 

■^n,k,t [^] -An,k,t [TZ+l] ^n,k,t [^] [T/] [/?Z+l] 

^ U 7i]i4 n> fe^[a n 7/]-^4n,/c,t [/^z+i] • 

Due to cm , a U 7; is a set of consecutive integers, so an application of cm yields 

LJ 7z]-^n,/c,£ [/^Z+l] ^ U 7/+l] -A-n,k,t [(^ U 7 z) D -|-1 ] - 

But (a U 7/) fl fii+i = a D (3i +1 since the sets & are pairwise disjoint. So, 

A- n ,k,t [TZ+l] — hJ 7Z+l]-^-n,A:,t D 7z]^n,/c,t D 

D / yi-\-l\An i k i t \p^ D 7/+l]. 

Now, given a set of consecutive integers a C (n) and any set /? C (n), write f 3 = 71U72 where 
7 i ; =U- = 1 A, 72:=U-1 T+iflii all & (i = 1,... ,1 + m) are separated, and Pi satishes (11.131) if 
and only if i < l. Then 

A-n t k,t [®] A n k i\P\ -4n,fc,t[o!]^n,fc,t[7l]-'4ra,fc,t[72] — C 7l]^-n,fc,t[® C 7l]-^n,A:,t [72] 

= A nAt [a U 71 U 7 2 ]A ntk> t[a fl 71] = A nAt [a U P\A nAt [a D 0\. 

In other words, A n ^,t satishes (11.111) if a C (n) is a set of consecutive integers and f3 C (n) 
is arbitrary. 

Finally, if 07, «2, P Q (n), the sets a* (z = 1, 2) are separated, (11.111) holds for 07 and /?, 
and «2 is a set of consecutive integers, then cm holds for a := 07 U «2 and /?: 

An,k,t A nk i [cri]v4 ni / c ^[cr2]v4 ri] fc i j [/3] 

fl C /?] A n k i [tti |7 

> A n)fcif [(o7 U /5) U a 2 ]A nifcif [(ai U /?) fl n p\ 

— A n) k,t[cv- U P]A n ^j[a.i n p]A n ^ it [a.2 n p] 

= A n ^t[a U P]A nAt [a fl P\. 

So, by induction on the number of ‘components’ of a, cm holds for any a, P C pi). 
Thus, an application of the Gantmacher-Krein-Carlson Theorem concludes the proof of the 
following. 

Proposition 2. The matrices A n ^ tt are GKK for any n, k & I N and f G (0,1). □ 
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1.3.3 A nt k,t are r-matrices 

Now check that have eigenvalue monotonicity for any n G IN and t G (0,1). If 

(n) I) a = U \ =1 OLi is the union of separated sets of consecutive integers, then 

o 

det (A nAt (a) - XI) = JJdet(H(cq) - XI) 

1=1 

since A n)k — XI is Hessenberg (the same observation earlier led to (11.101) 1. Since A U:kjt — XI 
is Toeplitz, the product in the right hand side equals J^ =1 D{# ai )(X). Hence, to prove 
eigenvalue monotonicity of A n k j it is enough to prove it for leading principal submatrices of 
A U)k ,t only, i.e., to show 

l(An+i,k,t) < l(A n ,k,t) Vn G IN. 

First recall that a P-matrix has no nonpositive real eigenvalues, since the coefficients of 
its characteristic polynomial strictly alternate in sign (Kellogg’s theorem gives a stronger 
statement, but one does not need it here). So, l(A ntkjt ) > 0 for all n G IN. Now note that 
D n {X) = (1 — A)" for n — 0,..., k+ 1. For n = k + 2, recall that a k+ \ = (—l) fc (l —f), so (11.81) 
yields D n {X) = (1 — A) n — (1 — t). Hence, 

l(Ak+2,k,t) = 1 — (1 — £)V( fc + 2 ) < l = l(A k+ i tk}t ) = • • • = l(Ai t k,t) 

and Djil^A^xt)) < 0 for all j < k + 2. 

On the other hand, if 

Dj(l(A nXt )) <0 Vj < n, (1.14) 

then m and the positivity of l(A n>k j) imply that 

fc+i 

D n +l(l(A n ,k,t)) — — (1 — t)KAn,k,t) Dn-j(l(A n ,k,t )) < 0 . 

3 = 1 

But P n+1 (0) = ff n ~ k B > 0, hence D n+1 changes its sign on the interval (0,/(H n fc t )), hence 
Z(H n+l fc t ) < l(A n)k j) and (11.141) holds for n + 1. Thus, by induction, the matrices A n ^ k)t have 
eigenvalue monotonicity. 

Proposition 3. A Ujkjt is a t- matrix for any n,fe G 1N and any f G (0,1). □ 

1.3.4 A 2k + 2 ,k,t is unstable for sufficiently large k and small t 

Now let B k := lim t _ > o+ A 2 k+ 2 ,k,t- The matrix B k is Toeplitz with first column 



2k terms 


and first row 

(1,0 1 _^ 1 0, (-l) fe , (-l) fc , 0_^0). 

k terms k —1 terms 
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Show that there exists K e IN such that, for all k > K , B k has an eigenvalue A with Re A < 0. 
As the eigenvalues depend continuously on the entries of the matrix, this will demonstrate 
that, for any k > K, there exists t G (0,1) such that the GKK r-matrix A 2k + 2 ,k,t das an 
eigenvalue with negative real part. 

The polynomial D 2k + 2 das a root with negative real part iff the polynomial defined by 

: = ( \ k ^X) k -^ = ^ + A) fc+3 — (k + 1)(1 + A) + k 
has a root with positive real part. Since 


fc+i 


*(A) = A[E(*i 3 

3 =0 V J 

it is, in turn, enough to show that r] k defined by 


A fc+3-;-i + 2 


k +3 


-fc(A) := ' j = 2A*« + g f + 3 I A l «-J 


has a root with positive real part. 

One is now in the position to apply the classical negative stability criterion of Hurwitz to 
the polynomial r] k . It is more efficient, however, to use the following condition necessary for 
nonpositive stability due to Ando. 

Theorem (Ando [T]). The Hurwitz matrix 


( 


a 3 

«5 

a 7 ■ 

•• \ 


ao 

a 2 

04 

Og • 



0 

CL\ 

03 

05 • 



0 

ao 

02 

04 • 


\ 





- ) 


dxd 


of a nonpositive stable polynomial f{x) =: Yl^=o a j x<i J °f degree d is totally nonnegative. 
The Hurwitz matrix for the polynomial r/ fc is 


H,:= 


(T) 

2 

0 
0 
0 

V 1 


e: 3 ) m m ra 
e: 3 ) m m m 
m m m m 

2 m m (W) 

o et 3 ) ( fc : 3 ) m 


/ 


(fc+2) X (fc+2) 


Compute the minor H k [ 2:5], taking out the factors ( k f 3 ), ( k ^ 3 ), ( fc g 3 ) from its second, third, 
and fourth columns respectively. This gives 

■mw 0 {3k3 ~ m2 - 210k ~ 318)(fc+4)2{k+5) (" 2 3 ) C V)( k t 3 ' 


H k { 2:5] = 


11 














Thus, Hk[ 2:5] < 0 for large enough k, precisely, for all k > 20. So, for k > 20, ip- has a zero 
with positive real part, therefore, D 2 k +2 has a zero with negative real part. 

This completes the proof of the following. 

Theorem 1 . The GKK t -matrices A 2 k+ 2 ,k,t are unstable for sufficiently large k and suffi¬ 
ciently small positive t. □ 


1.3.5 Numerics 


To illustrate the result, consider the matrix ^ 44 , 21 , 1 / 2 ) he., the Toeplitz matrix whose first 
column is 

(i,i,(V_o) T 

42 terms 

and first row is 

(1, 0 i _^0, -1/2, —1/2 2 ,1/2 3 , —1/2 4 ,..., -1/2 22 ) 

21 terms 

and the limit matrix B 21 , with the same first column as ^ 44 , 21 , 1/2 and first row equal to 




21 terms 20 terms 

According to MATLAB, the two eigenvalues with minimal real part of the first matrix are 
-2.809929189497896 • KT 2 ± 3.275076252367531 • l(Th; 
those of the second are 

-3.420708309454068 • 1(T 2 ± 3.400425852703498 • KTh. 


Further MATLAB calculations suggest that the matrices ^ 28 , 13 ,* are already unstable for 
sufficiently small t. 


1.4 More on the counterexample matrices 

1.4.1 The sign pattern of A~\ t 

The matrices A n ^,t turn out to be quite ‘close’ to being sign-symmetric matrices for n < 2k + 2. 
Namely, all their minors of order n — 1 are nonnegative. Since det A n ^ jt > 0 whenever t > 0, 
this property can be restated as follows. 

Proposition 4. A~ kt is checkerboard for n <2 k + 2. 

The proof will make use of the Gohberg-Semencul formula for inverting a Toeplitz matrix. 
Theorem (Gohberg and Semencul [I201J or, e.g., [221, p.21]). If the Toeplitz matrix T 
of order n is invertible, then 



/ Xo 

0 • 

■■ 0 \ 


/ yo y -1 • 

y~n +1 ^ 

T 

Xi 

Xo 

■■ 0 


0 yo ■ 

y~n +2 


\ x n _i x n _ 2 ■ ■ ■ x 0 ) \ 0 0 ■• ■ y 0 ) 
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/ 0 0 

V—n +1 o 

y~n +2 y~n +1 

V y-i y-2 


o 

o 

o 


°\ 
o 

0 


/ 0 
0 


1 

0 


%n—2 
%n—1 


^2 


\ 


.(1.15) 


y~n +1 




o 

o 


o 

o 


•£n—1 

0 


where x :=(xo,..., x n _i) T is £/ie solution to Tx = e\ and y :=(y_ n+1 ,..., j/ 0 ) T is the solution 
to Ty = e n (recall that e 3 denotes the j-th unit vector). 

Proof of Proposition f. There is nothing to prove if n < k + 1. The proof for n > k + 2 
is by induction on n. To avoid confusion, let us make the notation more precise. The goal 
is to inductively solve equations A n ^ it X n = e{ n) and A nt 
Let s := 1/t. Prove by induction that 


Xl 


yi 


( (—1 ) l s l+1 if l < n — k — 2 , 

| (— l) i s ri_fc ” 1 if l > n — k — 1, 

1 s if l — 0 

< (—l) z+1 (l — s) if — 1 > l > —k — 1 

lo if l < - 2. 


0,..., n — 1, 

l = —n + 1,..., 0 


(1.16) 

(1.17) 


for n > k + 2. Indeed, formulas (11.161) and (11.171) hold for n = k + 2 (which can be checked 
by direct calculation). To justify the inductive step, one needs to show that 


1 

pfa- 1 ) 


where B n ^ t 


B n ,k,t \ f S \ _ (n) ( 1 \ f 0 \ 

a,- w y v - sX n- t ) ~ ei • ! e'”- 11 / v / 

:= A n) fc )t (l, 2:n). Expanding (11.41) to order 2k + 2, one verifies that 



_ f (—l) fc (l — t) if j = k + 1 

aj ~ \ (-l)H j - k - 2 (l -t) 2 if k + 1 < j < 2k + 2 - 

So, 

n—k—1 n—k—1 

B n , k , t X n _! = ^ a i( -l) fe+i -V- fe - 2 = s ra - fe - 2 [ (1 - A) - ^ f- 2 (l-A) 2 ] = 

2=1 2=2 

1 4-71 —k —2 

s n - k ~ 2 [( 1 - A) - (1 - A) 2 ] = s n - fc - 2 (l - A )t n ~ k ~ 2 = (1 - A), 

(1 - A) 

so s(l - B n<ktt X n _i) = 1. Since e^ _1) = —l by the inductive hypothesis, this 

justifies the transition from X n _\ to X n . Now, 


Bn,k,A^n— 1 


7i—k—2 

ca(— l) n k *(1 — s) + a n _fc_is 

2=1 

7i—k—2 

(-l) n “ 1 [(l - s)((l - A) - (! - A)¥- 2 ) + t n ~ k ~\ 1 - A) 2 s] 

i=2 
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= (-ir 1 

= (-ir 1 


t -1 
t 

t - 1 
t 


-| j n—k—6 

(( i _ A )_( i _ A ) 2 ___) 

(1 - X)t n ~ k ~ 3 + t n ~ k ~ 3 ( 1 - A) 2 


+ t 


n—k—3 


l 

t 


= o 


(l-A) 2 


1 

t 


and, since 
This implies 


by the inductive hypothesis, the inductive step is completed. 


\xiy~m\ > \x n - m yi- n \ whenever m > 1, n-m> l. (1.18) 

Indeed, if l — n < —k — 2, the right hand side of (fT"T%l) is zero. But if Z > n — k — 1, then 
\xi\ = \x n _ m \ and \y_ m \ = \y t _ n \ (observe that l - n ^ 0). 

Since the (i,j) element of the right hand side of (11.151) has the form 

^ ^ %iVi—j—l ^ ^ %lUi—j—l 

l=i—min{i,j} l=n—j 


and the terms in both sums have sign (—1)* 3 , the subtraction gives 


f x 0 yi-j + (-1)* 3 r iij if i < j 
\ Xi-jyo + (-1 ) l ~ 3 r itj if i > j 


for some > 0. 


So, in either case, if the (i,j) element of A n \ t is nonzero, then its sign is (—1)* 3 ■ 


□ 


1.4.2 The spectrum of ^4oo,fe,t 


It is well known (e.g., m p.21]) that the spectrum of an infinite Toeplitz matrix T = ( r(i—j )) 
with 

OO 

J2\r(j)\<oo (1.19) 

j=-oo 


is the union of the curve 


C(T) = {S(T)(s) :s6 <D, |s| = 1} 


( 1 . 20 ) 


and those points in C \C(T ) whose winding number with respect to C{T) is nonzero. 
Proposition 5. The spectrum a(Aoo^j) is the union of 


C(^oo,fc,t) 


1 + ts 




:s6 C, 


1 } 


and the set of points in C whose winding number with respect to C(Ax>,fc,t) is nonzero. In 
particular, C(^4oo,fc,t) contains a negative real point. 

Proof. The value of the symbol d>(Aloo,fc,t) (given by (11.41) 1 at the point s = — 1 is 

-{l-t)/{k + 2-(k + l)t). (1.21) 
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□ 

Next, let us try to determine the limit set of eigenvalues of the finite sections of 

(the matrices A n ^,t) as n tends to infinity, ft is known (see, e.g., [TJ ) that £(T) C a{T) 
for any Toeplitz matrix satisfying ADD - Moreover, if S{T) is a rational function, there is a 
characterization of the set £(T) due to K. M. Day. 

Theorem (Day |12| h Let T : = (r(f — j))°° =0 be a Toeplitz matrix satisfying M.llA) with 
symbol S(T) that coincides with the expansion of the rational function -AF in the annulus 
{s G C : i?i < |s| < R 2 }, G (H) being a polynomial all of whose roots lie in the set 
|s| < R\ (\s\ > R 2 ), F a polynomial having no common factors with GH, and p:=degG. 
Let T n :=(r(f — j))^ =0 . Then the set 

CiT) :={A : A = limA m , A m € u(Tj m )} (1-22) 


coincides with 

(A € C : |r p (A)| = |r p+1 (A)|}, 

where rj( A) are the roots of the polynomial R := F — XGH listed in the order of their absolute 
values |r 1 (A)| < |r 2 (A)| < • • ■. 

For the problem in hand, F(s) = 1 + ts, G(s) = s, H(s ) = 1 + (1 — t) Xp=i ( — ■ s ) j 1 V = \ 
and R(s) = 1 +ts — As(l + (1 — t) Xp=i ( — S Y)- ^ seems hopeless to seek an explicit formula 
for £(^4oo,fc,t) for arbitrary k and t. Even checking whether a particular point belongs to 
£(Ax>,fc,t) turns out to be rather nontrivial. 

However, one can easily make the following (negative) observation. 

Proposition 6. The point \1.21\) does not belong to for any k 6 IN and t e (0,1). 

Proof. Let A equal (11.211) . If s < 1, then 

1^)1 >l-t, (1-23) 


while, by the triangle inequality, 


\XG(s)H(s)\ < |A|(1 + (1 - t)(k + 1)) = 1 -t, 


(1.24) 


so R = 0 iff inequalities (TT23P and dni both become equalities iff s = — 1. On the other 
hand, -^R(-l) 7 ^ 0, so s = — 1 is not a double root of R. Thus, |r*!(A)| < |r 2 (A)| and 
A ^ £(^4oo,fc,t)- O 

Nevertheless, some eigenvalues of A n ^t form sequences approaching points on the negative 
real axis as n tends to infinity. This can be shown using the following result of M. Biernacki. 
Theorem (Biernacki PJ as cited in [ill0 ]). Let p, q G IN be relatively prime and let 
f(s) := s~ p + s q — A. Then the set 

(A : |r„(A)| = |r p+ i(A)|} 

is the star-shaped curve 

S :={A = er} 

where e is any (p + q)th root of unity and 

0 < r < (p + g)p-P/(P+'?) g - 9 /(P+'?)_ 
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Here is the more precise formulation of the above claim. 

Proposition 7. For any k G IN, there exists t(k ) G (0,1) such that the set C(Aoo,k,t) contains 
a segment on the negative real axis for any t G (Q,t(k)). 

Proof. lit — 0, R(s ) :=(1 + s)R(s) = 1 + s — As(l — (—s) fc+2 ). For A ^ 0, the equation 
R(s) = 0 is equivalent to the equation 

(-s)" 1 + (-s) fc+2 - (J. = 0 where n := ,_^ 1/( fc +3) - (1-25) 

Taking e := 1 and applying Biernacki’s theorem amounts to solving the equation 

1 - A - (-A) 1/(fe+3) r = 0 (1.26) 

for some r G (0, ry) where ry : = (fc + 3)(/c + 2) _1 ^ fc+3 k Since ry > 2 for all k > 0, the left hand 
side of dnnp changes sign between (A =) — 1 and A = 0 for all r G (2,rfc), therefore, there 
exists a segment A:=(X min , \ max ) C (—1,0) such that |ri(A)| = |r 2 (A)| for the polynomial R 
whenever A G A. 

Without loss of generality, one can assume that R has no double roots when A G A, since 
that assumption can rule out only a finite number of A’s. So, one can assume that such 
points are already excluded from A. 

If, for some A G A, R had two distinct roots of the form s and e|s|, with e:=±l, that 
would imply that one of the triangle inequalities 

|(A — l)s + A(—s) fc+3 | < |A — 11 • |s| + |A| ■ |s| fc+3 
|(A — l)s + A(—s) fc+3 | > ||A — 1| • |s| - |A| ■ \s\ k+3 \ 

becomes equality, which is possible only if s k+ 2 is real, hence, by the condition R(s) = 0, 
s itself must be real. But the assumption that R has two roots s, — s G 1R leads to the 
contradictions 2 = 0 if k is even and (A — l)s = 0, hence A = 1, if k is odd. 

So, R has no roots of the form s, ±|s| for A G A. 

This implies, first of all, that none of the roots of R with minimum absolute value equals 
— 1 when A G A, so that r’i(A), r 2 (A) are also roots of R with minimal absolute value. More¬ 
over, the set of roots with the smallest absolute value consists of (possibly, several) distinct 
non-real conjugate pairs Sj, ~s]. Since the roots of an algebraic equation are continuous func¬ 
tions of the coefficients, one of those pairs must stay a pair of complex conjugate roots with 
smallest absolute value as t runs from 0 to t(k ) for some sufficiently small value t(k). □ 

To visualize the sets a(A 00: k t t) and C(A 00 ^,t), here are four figures, two for k:=3,t :=0.2 
and two for for k := 10, t := 0.4, drawn by MATLAB. 
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Figure 1. The curve C(^4 oo,m) (black) and the sets cr(A 50)fc;t ) (green), 
< 7(^100 ,k,t) (cyan), a(A 200 ,k,t) (blue), a(A 400i k,t) (magenta), 
cr(A 8 oo,k,t) (red). Here k = 3, t = 0.2. 



Figure 2. Figure 1 zoomed in the part of ^{A^^ 
around the negative real axis. 
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5 _i_i_i_i_i_ i _ i _ i _ i _ i 

-1.5 -1 -0.5 0 0.5 1 1.5 2 2.5 3 3.5 

Figure 3. The curve CO^oo,k,t) (black) and the sets a(A 50 ^ t ) (green), 
cr(Aioo ,k,t) (cyan), a(A 20 o,k,t) (blue), a(A 400 ,k,t ) (magenta), 
cr(A 8 oo,k,t) (red). Here k = 10, t = 0.4. 



Figure 4. Figure 3 zoomed in the part of cr(A 00j fc >t ) 
around the negative real axis. 
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1.5 Open problems 

The following questions appear to deserve further investigation in connection with the GKK 
T-matrix problem. 

1. Can the matrices be approximated by r-matrices that are strict GKK, i.e., P- 

matrices satisfying 

A[a, (3\A[f3, a] > 0 Va, (3 e (n), U f3 — 1? 

The matrices A Uj k,t themselves are not strict GKK. If the answer is no, 

la. Are strict GKK matrices positive stable? Are strict GKK and r-matrices positive 
stable? 

2. Given a, (3 C (n) with = #/3, call the number — #(a fl (3 ) the dispersal of 
the minor A[a,(3\. The counterexample from Section 1.3 shows it is not sufficient for 
stability of a P-matrix A that the inequalities 

A[a , 0\A[/3, a] > 0 

hold for all minors of dispersal < d 1. Carlson’s theorem asserts that the value d — n 
is sufficient for stability. What minimal value of the parameter d would guarantee 
stability? In particular, does that value depend on n? 

Also, here are two less directly related questions, which arose in the construction of the 
counterexample. 

3. Given n G IN and positive numbers ( Pa)n^aC(n) satisfying the generalized Hadamard- 
Fisher inequality (11.111) . when is there a matrix A such that A[a] = p a for all a? 

4. For a matrix A, let 

b • 

b f = c i : = 7FP J = 0,..., n. 

#a=j \j) 

When is it true that 

Cj > Cj-ic j+1 , j = 1,..., n - 1? (1.27) 

These inequalities are known for diagonal matrices with positive diagonal elements 
(e-g., [2D p.51]) and go back to Newton. Since the numbers Cj are invariant under 
similarity, Newton’s inequalities (f rm also hold for all diagonalizable matrices with 
positive real eigenvalues. Do the GKK matrices satisfy G2ZD? 
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Chapter 2 


Inverses of special matrices 


2.1 Bounded invertibility problem 


The most general setup of the problem of this chapter is the following. 
Let A be a collection of matrices such that 

inf min{|z| : z G cr(A)} > 0, sup ||A|| < oo 
^e-4 Ae.4 

for some norm || • ||. What conditions on (A,) imply 


( 2 . 1 ) 


sup || A 1 || < oo? 
AeA 


( 2 . 2 ) 


In the easy case when the order of matrices is bounded above, the conclusion o holds 
without any additional hypothesis. 

Proposition 8. Let A be a collection of complex matrices satisfying K2.1\) for some norm 
|| • || and such that 


sup order(A) < oo. 
AeA 


(2.3) 


Then \2.2f) holds. 

Proof. Without loss of generality, one can assume that all matrices A G A belong to 
C nxn where n := sup Ag ^ order(A). Indeed, just replace each A by 


A. diag(A, /( n -orderpi)) + )- 


Then (12.11) holds for the new collection A. Also, A satisfies (12.21) iff A satisfies (12.21) . Next, 
since all norms on C nx?l are equivalent, one can assume that || ■ || is an operator norm 
subordinate to a norm on C ri (also denoted by || ■ ||). 

So, suppose A C C nx?l and (12.21) is violated. Then, by the Banach-Steinhaus theorem, 
there exists v G C n and a sequence (A,) such that 



oo 


But the sequence (A,-) is totally bounded in <D nxn , hence contains a Cauchy subsequence. 
Without loss, it is (A,-) itself. The limit A linij^oc A 3 is invertible since 


min{|z| : z G cr(A)} = lim min{|^| : 0 G cr(Aj)}, 
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hence lim ||A- 1 t;|| = ||A 1 n|| < oo. This contradiction shows that Q) holds for any operator 
norm, hence for any norm on C nxm . □ 

Next, let us consider the case when the order of matrices Aj is not bounded above and 
the norm in question is the oo-norm. This question was posed by K. West [31] for positive 
definite Hermitian matrices in connection with a problem from econometrics. As is shown 
in Section 2.3, the answer in that case is no. However, under certain additional hypotheses 
the answer is yes, as is shown next. 

2.2 Boundedly invertible collections of Hermitian matrices 

One of the possible restrictions on the collection (A,) that ensures that (12.21) holds is (uni¬ 
form) bandedness of matrices Aj, as follows directly from a theorem of S. Dernko. To recall, 
a matrix A = ( a{i,j )) is called banded with band width w if 

a(i,j ) = 0 whenever \i — j\ >w. 

Theorem (Demko [Qj). Let A £ C nxri be banded with band width w and satisfy conditions 
||A|| P < 1, ||A _1 || P < /i' 1 for some 1 < p < oo and p > 0. Then, with A -1 =:(a(i, j))’j =1 , 
there are numbers K and r £ (0,1) depending only on /a and u ; such that 

|«(bi)l < Vi,j = 1 ,...,n. 

In particular, for any 1 < q < oo, ||A _1 || 9 < C where the bound C depends only on /i, and 

w. 

Since the smallest eigenvalue of a Hermitian matrix is the reciprocal of the 2-norm of its 
inverse, the hypothesis of Demko’s theorem is satisfied for p:=2 whenever the collection A 
satisfies m , hence the collection A 1 is bounded in any g-norm (1 < q < cx)), in particular, 
in the oo-norm. 

A different restriction that ensures boundedness of (Aj 1 ), provided that A/s are Hermi¬ 
tian and satisfy m, is the oscillatory property. Recall that a matrix is totally positive if all 
its minors are positive. A totally nonnegative matrix A is called oscillatory if A 1 is totally 
positive for some l £ IN. It is well known (see, e.g., m p. 123]) that cr(A) consists of n 
distinct positive real numbers Ai < • ■ • < A n and that the fcth eigenvector Vk (Avk = XkVk) 
(unique up to a scalar multiple) has no zero entries and precisely n — k sign changes if A is 
oscillatory. 

Theorem 2. Let A 6 (p nxn y e an oscillatory Hermitian matrix with smallest eigenvalue 
A Th pt) 

/x mm • - 1 - 1 

IIA-'lloo < MUrR. (2.4) 

^min 

The proof will make use of the following lemma due to C. de Boor. 

Lemma (de Boor |BJ). Let A £ (j2 n><n y e a nonsingular totally nonnegative matrix. If, for 
some x, y £ C n , 

Ax = y, signrr(f) = signy(f) = (—1)* _1 , 
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then 


HA - 1 


< 



mirij \y(i)\' 


Proof of Theorem 2. Since A is Hermitian, its eigenvector v corresponding to the eigen¬ 
values A min can be uniquely determined from the minimization problem 

v*Av —> min 


once one of the entries of v is fixed. Let r>(l) = 1, so that v = [1 v\ T . Since A is real, all the 
entries of v are necessarily real. Let A be partitioned conformably to v: 


( a( 1,1) A(l,2:n)\ 

\ H(l,2:n) T A{2:n) )' 


Then v T Av = a(l, 1) + 2^4(1, 2 :n)v + v 1 A(2\n)v achieves its minimum at 

v := —A(2:n)^ 1 A(l, 2 :n) T . 

By the eigenvalue interlacing property of A, mincr(H(2:n)) > A m ; n . Since ||H(1, 2:n) T || < 
\\A\h, this yields 

INloo < IN 2 < P(2:n)- 1 || • ||H(l,2:n) T || 2 < 

The same argument can be applied to the case when any other entry of v is set to be 1. 
Since the eigenvector v is unique up to multiplication by a scalar, this means 

t~~ jU < whenever Av = X min v. 

nun i |u(?)| X min 

Now apply de Boor’s lemma with x : = v, y := X min v and obtain (ED . □ 


2.3 On shifted Hilbert matrices and their companions 

However, the Hermitian property alone is not sufficient for the implication , as 

is demonstrated by the two examples below. For those counterexamples, one needs several 

additional notions. The matrix H n := ( 1 , ) is known as the Hilbert matrix. It has been 

I'+J-VyU 

a subject of extensive studies and has served as an example of many unusual phenomena in 
operator theory. M.-D. Choi in P used what he called the companion of the Hilbert matrix , 

viz. the matrix C n := .It turns out that the matrix C n belongs to the special 

class of ultra- m etric matrices 1 introduced by Nabben and Varga as a generalization of the 
notion of strictly ultrametric matrices attributed by Martinez, Michon, and San Martin E2 
to C. Dcllacherie Pl 2 - _ 

1 This matrix is also from the class of single-pair or (‘one-pair’) matrices due to Gantmacher and Krein FT71 p. 113]. 
2 Ultrametricity was first introduced in connection with p-adic number theory. Recall that a distance d on a space X is said 
to be ultrametric if it satisfies the inequality 

d(x , y) < max{d(r, z), d(z, y )} \/x, y, z E X. 
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A matrix A—:(a(i,j))fj =1 is ultrametric if 
A = A*, A>0 

a(i,j) > min{a(i, k), a(k,j)} V i,j,k G (n) 
and 

a(i,i) > ma x{a(i,k) : k G (n)\{i}} Vi G (n). (2.5) 

If the inequality EH3> is strict for all i G (n), A is called strictly ultrametric. 

Finally, before stating the result of Martinez, Michon, and San Martin, recall that a 
matrix A—:(a(i,j)) G <D nxn is row (column) diagonally dominant if 

a(M) > Hhj)\ e ( n ) (2.6) 

(a(i,i) > E iG (n>\{i} l a CM)l Vi G (n)). (2.7) 

If all the inequalities (12.61) (the inequalities (12.71) 1 are strict, A is called strictly row (column) 
diagonally dominant. 

Theorem (Martinez, Michon, and San Martin [2?]). The inverse of a strictly ultra¬ 
metric matrix is a symmetric strictly diagonally dominant M-matrix. 

Now one can construct the following counterexample to the implication (EH) =► E2D- 
Proposition 9. Let a > 0 and let A n : = ( al n + C n ) _1 . Then the collection {A n )neN of 
Hermitian positive definite matrices satisfies (TO) but does not satisfy (EH) . 

Proof. Subtracting of the jth column from the j — 1st column of the matrix C n , for j = 
2,... ,7i, one verihes that det C n > 0 for all n G IN. Hence C n is a Hermitian positive def ini te 
matrix. So, min a(aI n +C n ) > a. By |TJ[ Problem V], ||C„|| 2 < 4, so ma xcr(aI n +C n ) < 4+a. 
Thus, <r(A n ) C [1 /(a + 4), 1/a] for any n G IN. 

The matrices C n are ultrametric, hence the matrices al n + C n are strictly ultrametric, 
so by the theorem of Martinez, Michon and San Martin, their inverses A n are diagonally 
dominant M- matrices. But any diagonal entry of A n is bounded above by 11 A n \\ 2 < 1/a, so 
the oo-norm of A n is bounded above by 2/a. On the other hand, 

P^lloo = II aln + C'nlloo ~ a + In n -KX). 

n —► oo 

So, the collection (A n ) satisfies (12.1|) but violates Mb □ 

Proposition 10. For large enough a > 0, the collection (A n :=al n + H n )~ l satisfies \2.1\) 
but does not satisfy 

Proof. Note that C n — H n > 0 and estimate ||C'„ — 77 n ||oo. 



By Dellacherie’s definition, a symmetric matrix A E C nXn is ultrametric if there exists an ultrametric distance d on (n) such 
that 

d(i, j) = d(i, k) iff a(i, j) = a(i, k). 
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Recall that 


P -1 || - 1-\\A-bI ||R-1]| whenever \\ A - B \\ ■ ll s_1 ll < 1 

for any operators A, B and operator norm || • ||. The proof of Proposition 9 demonstrated 
that || (al n + Cn)^ 1 ||oo < 2/a. Hence, if a > 4, then 

\\C n — H n \\oo\\{0iIn -\-C n ) 1 || 00 5 ; — < 1 ) 

a 

hence \\(al n + H n ) 1 Thus, the collection (A n = (a/ n + H n ) _1 ) satishes (12.111 but 

violates (JO- ° □ 


2.4 Inverses of nonnegative Hermitian Toeplitz matrices 


The last example in the same spirit deals with Hermitian Toeplitz matrices. 

Proposition 11. Let T k be the infinite (upper triangular) Toeplitz matrix with symbol 
S{Tk)(s): = 1 + s + cs k , where c is any complex number with |c| > 2 (chosen to be positive if 
the matrices A n must be nonnegative). Set A k : = T k Tf and let A njk be the leading principal 
submatrix of A k of order n. Then the collection (A ntk ) of positive definite Hermitian Toeplitz 
matrices satisfies (HH but violates \2.ty) . 

Proof. By the spectral theory of Toeplitz matrices, which was already discussed in 
Section 1.5, cr(A k ) equals the set of values of its symbol on the unit circle |s| = 1. Notice 
that | S{T k )(s )| > |cs fc | — |1 + s\ > \c\ — 2 > 0 whenever |s| = 1, so 

inf miner (A k ) > 0. 

fceiN 

On the other hand, the matrices A k have at most 9 nonzero diagonals with the absolute 
value of each nonzero term at most 2, so 


sup PfcHoo < 18. 
fce n 


Finally, Af 1 


T k l T k 1 and the symbol S(T k of T k 1 is 


1 

S(T k ){s) 


l-s + s 2 -+ (-l) fc_1 s fc - 1 H-. 


So, the {k — 1) x {k — 1) leading principal submatrix of A k l has the form 


( 1 

-1 

1 ••• 

(-l) fc \ 

-1 

2 

-2 ••• 

(-1)^2 

1 

-2 

3 ••• 

(—l) fc_2 3 

V (-i) fc 

(—l) fc “ 1 2 

(—l) fc_2 3 ••• 

k — 1 / 


hence sup fee]N ||4 fc 1 ||i = oo. 

Since the matrices A k are Hermitian, the limit set of the eigenvalues of A n ^ k (as n tends 
to infinity) coincides with a(A k ) (see, e.g., 0). Since the (elementwise) limit of A~\ is A k l 
by da P-74], the collection (A n k ) satishes (12.11) but not (J22J. ’ □ 
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2.5 Least-squares spline projection matrices 

An interesting problem of the same type arises in spline theory. 

Bounding the oo-norm of the (L 2 -)orthogonal projector onto splines leads to matrices of 
the specific form 

k f 

A n ,k,t(h j)' = - -— / BikBjk. ( 2 . 8 ) 

ti-\-k J 

Here n, k G IN, t is a nondegenerate knot sequence with n + k knots, B ik is the i-th B- 
spline of order k for the sequence t, S k ,t- = span{T> lfc ,.. ., B nk }, and Lf is the least squares 
approximation to / G L 00 \ti 1 t n+k \ by elements of S k)t • C. de Boor [5] showed that there 
exists a positive constant C k such that 

C k \\A~i t \U < ||L|U < IKltlU 

and conjectured that, for k fixed, 

SU PPnJ,tll°o < oo. (2.9) 

n,t 

This conjecture was recently proved by A. Shadrin jnj using sophisticated tools from spline 
theory to construct vectors x and y appearing in de Boor’s lemma with the min \y(i)\ and 
max |x(f)| depending only on k but not on the knot sequence t or n and conclude, using the 
lemma, that m holds. 

The matrices A n , k t are known to be oscillatory and diagonally similar to Hermitian ma¬ 
trices, i.e., 

An,k, t D n , k ,tHn,k,tD nk ^ 

where D n k t are diagonal with positive diagonal entries. (For sure, the condition number 
of the matrices D n k t is not bounded above.) Moreover, t are banded with band width 
k. Finally, the smallest eigenvalue of H Hyk}t is known to be bounded away from zero inde¬ 
pendently of t and n, so the collection (A n>k}t ) satisfies (12.11) . (All those facts can be found 
in, e.g., US p.401-406].) However, the above conditions are not sufficient for (12.21) (that 
is, (12.91) 1 to hold. In particular, an application of Theorem 2 to the eigenvector v m i n of H rhk t 
corresponding to its smallest eigenvalue demonstrates that the eigenvector D n>kjt v of A Uik ,t 
corresponding to the same eigenvalue cannot be used to prove (EH. 

Hence the following (somewhat vaguely formulated) problem: 

Problem. Single out an additional property of the matrices A n ^ k)t to obtain a simple matrix 
theoretic proof of \2. y\) - 


25 



Index 


Ando’s theorem, 13, 41 

Banach-Steinhaus’ theorem, 30 
banded matrix, 30, 38 
Biernacki’s theorem, 20, 41 
de Boor’s conjecture, 38, 41, 44 
de Boor’s lemma, 32, 33, 38, 41 

Carlson’s conjecture, 3 
Carlson’s theorem, 4, 27, 42 
Cauchy-Binet formula, 4 
checkerboard matrix, 15 
condition number, vi, 38 

Day’s theorem, 18, 19, 42 
Demko’s theorem, 31, 42 
diagonal matrix, 4, 28, 38 
diagonal scaling, 4 
diagonally dominant, 34, 35 
strictly, 34 
dispersal, 27 

eigenvalue interlacing, 2, 32, 43 
eigenvalue monotonicity, 2, 3, 11, 12, 42 
Engel-Schneider’s conjecture, 3, 42 

Gantmacher-Krein-Carlson’s theorem, 8, 10, 
41, 42 

GKK matrix, ii, 1-5, 8, 10, 12, 14, 27, 28, 
43 

strict, 27 

Green’s matrix, 33 

Hadamard-Fisher inequality, 8, 27, 42 
Hermitian matrix, ii, iii, 2, 3, 30-34, 36-38 
positive definite, ii, 2, 3, 30, 34, 36 
Hershkowitz’ conjecture, 4, 43 
Hessenberg matrix, 5, 8, 11 


Hilbert matrix, 33 
companion of, 33 
Hurwitz matrix, 13 

Kellogg’s theorem, 4, 11, 43 

Laplace expansion, 5, 7 
Laurent series, 6, 42, 44 

M-matrix, 1-3, 34, 35, 44 
Martfnez-Michon-San Martin’s theorem, 34, 
35, 43 

matrix norm, ii, v 

Newton’s inequalities, 28 

one-pair matrix, 33 
operator norm, v, 30, 35 
oscillatory matrix, ii, 31, 38 

P-matrix, 1, 2, 4, 8, 11, 27 
p-norm, v 

2-norm, 31 

oo-norm, ii, 30, 31, 35, 37 
Perron-Frobenius theory, 2 

sign-symmetric, 1, 2, 4, 15 
weakly, 1, 41 
spectral radius, v, 1 
spectrum, ii, v, 3, 18 

of a Toeplitz matrix, 18 
stability, ii, 3, 4, 13, 27, 42, 43 

T-matrix, ii, 2, 4, 5, 11, 12, 14, 27, 43, 44 
Taylor series, 6 

Toeplitz matrix, ii, iii, vi, 5, 6, 11, 12, 14, 
15, 18, 19, 36, 41-44 
symbol of, 6, 19, 36 


26 



totally nonnegative matrix, ii, 1-3, 13, 31, 
32 

totally positive matrix, 31, 41, 43 

ultrametric, 33, 35, 43, 44 
strictly, 33-35 
unit vector, iv, 16 

Varga’s conjecture, 3, 44 

West’s problem, 30 


27 



Bibliography 


[1] T. Ando, Totally positive matrices, Linear Algebra Appl. 90 (1987) 165-219. 

[2] A. Berman, R. Plemmons, Nonnegative matrices in the mathematical sciences. SIAM, 
Philadelphia, 1994. 

[3] R. Bhatia, Matrix analysis. Springer-Verlag, New York, 1997. 

[4] M. Biernacki, Snr les equations algebriques contenant des parametres arbitraires, Bull. 
Acad. Polon. Sci. (1927) 541-685. 

[5] C. de Boor, The quasi-interpolant as a tool in elementary polynomial spline theory, 
in “Approximation Theory” (C. G. Lorentz, ed.), pp. 269-276, Academic Press, New 
York, 1973. 

[6] C. de Boor, On a max-norm bound for the least-squares spline approximant, in “Ap¬ 
proximation and Function Spaces,” Proceedings, International Conference, Gdansk, 
1979, pp. 163-175, New York, 1981. 

[7] A. Bottcher, S. M. Grudsky, B. Silbermann, Norms of inverses, spectra, and pseu¬ 
dospectra of large truncated Wiener-Hopf operators and Toeplitz matrices, New York 
J. Math. 3 (1997), 1-31 (electronic). 

[8] D. Carlson, weakly sign-symmetric matrices and some determinantal inequalities, Col- 
loq. Math. 17 (1967) 123-129. 

[9] D. Carlson, A class of positive stable matrices, J. Res. Nat. Bur. Standards Sect. B 78 
(1974) 1-2. 

[10] A. Cauchy, Sur l’equation a l’aide de laquelle on determine les inegalites seculaires de 
mouvement de planetes, Oeuv. comp. 9(2) (1829), 174-195. 

[11] M.-D. Choi, Tricks or treats with the Hilbert matrix, Arner. Math. Monthly 90 (1983), 
no.5, 301-312. 

[12] K. M. Day, Toeplitz matrices generated by the Laurent series expansion of an arbitrary 
rational function, Trans. Amer. Math. Soc. 206 (1973) 224-245. 

[13] C. Dcllacherie, Private communication (1985). 

[14] S. Demko, Inverses of band matrices and local convergence of spline projections, SIAM 
J. Numer. Anal. 14 (1977), no.4, 616-619. 


28 



[15] R. DeVore and G. Lorentz, Constructive approximation. Springer-Verlag, Berlin, 1993. 

[16] G. M. Engel and H. Schneider, The Hadamard-Fisher inequality for a class of matrices 
defined by eigenvalue monotonicity, Linear and Multilinear Algebra 4 (1976) 155-176. 

[17] F.R. Gantmacher, M.G. Krein, Oscillation matrices and kernels and small vibrations 
of mechanical systems. Gostechizdat, 1950. 

[18] I. C. Gohberg and I. A. Fel’dman, Convolution equations and projection methods for 
their solution. Translations of Mathematical Monographs, vol. 41, AMS, Providence, 
1974. 

[19] S. Friedland, Weak interlacing property of totally positive matrices, Linear Algebra 
Appl. 71 (1985) 95-100. 

[20] I. Gohberg and A. A. Semencul, On inversion of finite-section Toeplitz matrices and 
their continuous analogues (in Russian), Matem. Issled., Kishinev, 7(2) (1972) 201-224. 

[21] G. Hardy, J. E. Littlewood, and G. Polya, Inequalities, Cambridge LIniv. Press, 
Cambridge-New York, 1952. 

[22] G. Heinig and K. Rost, Algebraic methods for Toeplitz-like matrices and operators, 
Birkhauser Verlag, Basel-Boston, Mass., 1984. 

[23] D. Hershkowitz, Recent directions in matrix stability, Linear Algebra Appl. 171 (1992) 
161-186. 

[24] D. Hershkowitz and A. Berman, Notes on u- and r-matrices, Linear Algebra Appl. 58 
(1984) 169-183. 

[25] O. Holtz, Not all GKK r-matrices are stable. Linear Algebra Appl. 291 (1999) 235-244. 

[26] R. B. Kellogg, On complex eigenvalues of M and P matrices, Numer. Math. 19 (1972) 
170-175. 

[27] S. Martinez, G. Michon, and J. San Martin, Inverses of ultrametric matrices are of 
Stieltjes type, SIAM J. Matrix Anal. Appl. 15 (1994) 98-106. 

[28] V. Mehrmann, On some conjectures on the spectra of r-matrices, Linear and Multilin¬ 
ear Algebra 16 (1984) 101-112. 

[29] R. Nabben and R. Varga, Generalized ultrametric matrices - a class of inverse M- 
matrices, Linear Alg. Appl. 220 (1995) 365-390. 

[30] P. Schmidt and F. Spitzer, The Toeplitz matrices of an arbitrary Laurent polynomial, 
Math. Scand. 8 (1960), 15-38. 

[31] A. Shadrin, The Loo-norm of the L 2 -spline-projector is bounded independently of a 
knot sequence. A proof of de Boor’s conjecture, manuscript (1999). 

[32] O. Taussky, Research problem, Bull. Amer. Math. Soc. 64 (1958) 124. 


29 



[33] R. Varga, Recent results in linear algebra and its applications (in Russian), in: Nu¬ 
merical Methods in Linear Algebra, Proceedings of the Third Seminar of Numerical 
Applied Mathematics, Akad. Nauk SSSR Sibirsk. Otdcl. Vychisl. Tsentr, Novosibirsk, 
1978, pp. 5-15. 

[34] K. West, Private communication (1999). 


30 



